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Abstract. We use numerical simulations to test a broad range of plausible observational strategies designed 
to measure the time delay between the images of gravitationally lensed quasars. Artificial quasar light curves 
are created along with Monte-Carlo simulations in order to determine the best temporal sampling to adopt 
when monitoring the photometric variations of systems with time delays between 5 and 120 days, i.e., always 
shorter than the visibility window across the year. Few and realistic assumptions are necessary on the quasar 
photometric variations (peak-to-peak amplitude and time-scale of the variations) and on the accuracy of the 
individual photometric points. The output of the simulations is the (statistical) relative error made on the time 
delay measurement, as a function of 1- the object visibility over the year, 2- the temporal sampling of the light 
curves and 3- the time delay. Also investigated is the effect of long term microlensing variations which must be 
below the 5% level (either intrinsically or by subtraction) if the goal is to measure time delays with an accuracy of 
1-2%. However, while microlensing increases the random error on the time delay, it does not significantly increase 
the systematic error, which is always a factor 5 to 10 smaller than the random error. Finally, it is shown that, 
when the time delay is comparable to the visibility window of the object, a logarithmic sampling can significantly 
improve the time delay determination. All results are presented in the form of compact plots to be used to optimize 
the observational strategy of future monitoring programs. 

Key words. Gravitational lensing: time delay, quasar, microlensing - Cosmology: cosmological parameters, Hubble 
constant. 


1. Measuring time delays 

Measuring time delays in gravitationally lensed quasars 
is difficult, but not as difficult as it first appeared in the 
late 80s when the first monitoring programs were started. 
Obtaining regular observing time on telescopes in good 
sites was (and is still) not easy and the small angular 
separations between the quasar images require to perform 
accurate photometry of blended objects, sometimes with 
several quasar images plus the lensing galaxy within the 
seeing disk. 

1.1. COSMOGRAIL 

The COSMOGRAIL project (COSmological MOnitoring 
of GRAvltational Lenses), started in April 2004, addresses 


both issues of carrying out photometry of faint blended 
sources and of obtaining well sampled light curves of 
lensed quasars. The project involves 5 telescopes: ( 1 ) the 
Swiss 1.2m Auler telescope located at La Silla, Chile, (2) 
the Swiss-Belgian 1.2m Mercator telescope, located in the 
Canaria islands (La Palma, Spain), ( 3 ) the 2m robotic 
telescope of the Liverpool University (UK), also located 
at La Palma, ( 4 ) the 1.5m telescope of Maidanak obser¬ 
vatory in Uzbekistan, and ( 5 ) the 2m Himalayan Chandra 
Telescope (HCT). 

All 5 telescopes, and others that will join the collab¬ 
oration, are used in order to follow the photometric vari¬ 
ations of most known gravitationally lensed quasars that 
are suitable for a determination of Hq . The sample of tar¬ 
gets is described further in Saha et al. (2005), as well as 
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the non-parametric models and predicted time delays for 
all of them. Our target accuracy on individual photomet¬ 
ric points is 0.01 mag, well within the reach of a l-2m class 
telescope and average seeing (1") in a good site. This accu¬ 
racy is reached even for the blended components of lensed 
quasars, thanks to image deconvolution algorithms such 
as the MCS algorithm (Magain et al. 1998). 

The results presented in the following were obtained to 
optimize the observations of the COSMOGRAIL project, 
which aims at measuring time delays, with an accuracy 
close to 1% within 2 years of observations (per lensed 
quasar). 

Although large amounts of time are available for 
COSMOGRAIL on each telescope, optimizing the time 
spent on each lensed quasar, depending on its luminos¬ 
ity and expected time delay, remains very important. The 
present paper aims at optimizing the temporal sampling 
to adopt in order to derive accurate time delays for as 
many lensed quasars as possible. 

The paper is organized in the following way. Section 2 
describes how we simulate the light curves of the quasar 
images. In Section 3, we present which parameters of the 
simulated light curves are varied and in which range they 
are varied. In Section 4, we explain how the time delays 
are extracted from the simulated light curves. The results 
of these simulations are discussed in Section 5. Since most 
lensed quasar light curves are probably affected by mi- 
crolensing events, it is important that our simulations in¬ 
clude such effects in order to evaluate their influence on 
the determination of the time delay. This is treated in 
Section 6. Finally, Section 7 investigates the effect of log¬ 
arithmic sampling on the light curves and shows how this 
irregular sampling can improve the time delay measure¬ 
ments when it is of the order of the visibility window 
of the object. Note that we consider here only the time 
delays measured from optical light curves. Radio obser¬ 
vations have characteristics that are completely different 
from the present simulations: noise properties, better spa¬ 
tial resolution, less sensitivity to microlensing. 
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Fig. 1. Expected quasar variations Sm in magnitudes as a 
function of the time interval St (see text). The curves are 
plotted for 2 different absolute magnitudes Mb and for 3 
quasar redshifts. 

wavelengths. Most other lensed quasars have time delays 
known with a precision of about 10%. 

The accuracy of the time delays depends critically on 
the temporal sampling, on the visibility of the object over 
the year, on the influence of microlensing, and on the good 
will of the quasar source to show photometric variations at 
all. Using numerical simulations on artificial quasar light 
curves, we try in the present work to define the optimal 
observational strategy to adopt in order to reach a desired 
accuracy on the time delay. We consider only the time de¬ 
lay between two quasar images. Our simulations remain 
applicable to multiple time delays in quads, but the er¬ 
rors on the photometric measurements of the 4 (or more) 
components must be uncorrelated. 


1.2. Which accuracy 7 

Not all lensed quasars are suited to an accurate determina¬ 
tion of Ho, first because not all of them have nice lens mod¬ 
els with little influence of degeneracies and, second, be¬ 
cause the error on the time delay propagates linearly into 
the error budget on Hg. While the latter is not the domi¬ 
nant component in the error budget it can (and should) be 
made almost negligible compared with the other sources 
of uncertainty. A precision of a few percent should be the 
goal of the photometric monitoring programs aimed at 
measuring time delays, if Hq is to be measured with an 
accuracy competitive with other methods. So far, very few 
time delays are known with very high accuracy. Among the 
best examples are the double Q 0957-1-561 (Colley et al. 
2003), measured in optical wavelengths, and the quadru¬ 
ple B 1608-1-656 (Fassnacht et al. 2002), measured in radio 


2. Simulated light curves 

The first step of the process is to generate artificial quasar 
light curves whose properties mimic quasar variations in 
a realistic way. A useful benchmark here is the analysis 
by Hook et al. (Hina, of the variability properties of a 
sample of some 300 quasars. 

They find that the variability 5m of optically selected 
quasars can be represented by a function of the form: 


\Sm\ = (0.155 -k 0.023 (Mb -k 25.7)) 



0.18 


where Mb is the absolute B magnitude, 2 ; is the redshift, 
and St is the time interval, in days, over which 5m ap¬ 
plies. We plot this relation in Fig. ^ for different lumi¬ 
nosities {Mb = —26, —28), and redshifts {z = 1, 2, and 
3), which are typical for lensed quasars. The curves in the 
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Fig. 2. Example of simulated light curve, for a 2-year long 
observation and a peak-to-peak amplitude A = 0.1 mag. 
The continuous light curve is shown as a solid line. It 
has been smoothed on a length scale of 30 days. The four 
samplings used in the simulations are shown (plus the log¬ 
arithmic sampling, see text), along with the error bars of 
0.01 mag. The figure is constructed for an object with a 
visibility of 8 consecutive months, hence the size of the gap 
in the center of the curves is 4 months. The curve plotted 
for the logarithmic sampling has the same number of data 
points as the curve for the 7-day sampling. 


figure show a variation time-scale of 10 to 100 days for a 
typical change of a few tenths of a magnitude. Realistic 
simulated light curves should show a time-scale for the 
variations, and a total “peak-to-peak” amplitude A = 5m 
of the variation that are in accord with these results. 

In order to mimic these variations we first consider 
a time series of N points spanning the total duration of 
the observation. For each of these points we define a sim¬ 
ulated magnitude in such a way that the artificial light 
curve follows a random walk with an arbitrary peak-to- 
peak amplitude. We next smooth the light curve with 
a Gaussian kernel that has a Full-Width-Half-Maximum 
(FWHM) of 30 days to ensure that the typical variation 
time-scale matches that of real quasars. Finally, the curve 
is renormalized so that its maximum variation is equal to 
a specified peak-to-peak amplitude A. 

A second light curve is then obtained, by applying a 
time shift > 0. The two sets of points {Ia, mA) for 
image A and (<b, tob) for image B obtained in this way 
are the final simulated light curves, both sampled with 
10 points per 24 hours. This sampling, very small com¬ 
pared with the sampling that will be adopted to carry out 
the actual observations, ensures that no interpolation is 
necessary when shifting curve B relative to curve A. The 


precision on the shift is, then, 0.1 day, 50 times smaller 
than the smallest time delay we wish to simulate. 

The curves are used to produce artificial observations, 
this time with a much sparser sampling. We define Nobs 
observing points at the observing dates tabs- For each of 
these dates we define the observed magnitude by selecting 
the closest value in time among the pairs (tA, mA) for 
image A and (ts, tob) for image B, resulting in noise free, 
sampled, artificial light curves. Finally, simulated photon 
noise is added to the data. This is achieved by adding to 
each observing point a normally distributed deviate of zero 
mean and variance Cohs • Thus one has defined a combined 
set of Nobs observations, (tobs, 'mA,obs) for image A and 
(tobs, 'mB,obs) for image B. Typical light curves are shown 
in Fig. El 

3. Parameter space 

In the simulations presented below, some parameters are 
imposed on us by technical limitations. This is the case 
of the maximum accuracy on the photometry of the in¬ 
dividual quasar images. We assume that a good goal is 
0.01 mag for a typical lensed quasar, or a signal-to-noise 
of 100 integrated over the quasar image. We have tested 
some cases where the points have larger error bars, and 
this led to the conclusion that the adopted 0.01 mag error 
is a requirement to meet in order to carry out the project 
successfully. Doubling the error bars also doubles the er¬ 
ror on the time-delay determination. Errors above 0.05 are 
likely to compromise the whole feasibility of the project. 
We also suppose that the algorithm used to carry out the 
photometry on the real data actually yields photon noise 
limited measurements. Second, we fix the total duration 
of the observations to two years, since one probably wants 
to measure Hg in a reasonable amount of telescope time. 

Other parameters cannot be fixed in advance. They 
define the parameter space we want to explore through 
the simulations, and include: 

1. The temporal sampling of the curves. We consider reg¬ 
ularly spaced sampling intervals of 3, 7, and 15 days 
(±30% due, e.g., to bad weather). Also, in some obser¬ 
vatories, large chunks of time are allocated rather than 
regularly spaced dates. To model this, as an example, 
we also carry out our experiment with a sampling of 
one observing point taken every other day during 15 
days, followed by 1 single point taken in the middle 
of the next month, and again one point every other 
day for 15 days, and so on. We refer to this type of 
sampling as “irregular sampling”. 

2. The visibility of the object. An equatorial object is seen 
no more than 5 months in a row in good conditions. 
A circumpolar object is by definition visible the entire 
year. We also choose an intermediate visibility of 8 
months. It should be noted that we do not allow for 
large losses of data points, e.g. non allocation of time to 
the project during a full semester, which would simply 
hamper even a rough estimate of the time delay. 
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3. The amplitude A of variation of the quasar. We choose 
three typical peak-to-peak variations oi A = 0.1, 0.2, 
0.3 mag over the two years of simulated observations, 
as suggested by Fig. Q] 

Determining the best combination of these parameters 
is the goal of the present work, for a broad range of time 
delays, from 5 to 120 days. For each time delay there are 
4 temporal samplings x 3 visibilities x 3 amplitudes = 36 
different possible combinations of parameters. 

4. Extracting the time delay 

Using the light curves constructed in the previous sec¬ 
tion, we now try to recover the time delay Atm chosen in 
the simulated data. Many cross-correlation techniques are 
available for this task, with a variety of technical subtleties 
dealing with unstable solutions, sparse sampling, and the 
effects of additional perturbations to the light curves (such 
as those caused by microlensing). 

The aim of the present experiment is to decide which 
observing strategy will assure us that the present typical 
10% error bar on optical time delays decreases below 2%, 
rather than testing the cross-correlation techniques them¬ 
selves. For this reason, without further discussion, we have 
adopted the cross-correlation method of Pelt et al. (Cnni, 
which is in wide use, and which combines robustness, sim¬ 
plicity, and low cost in terms of computing time. No other 
correlation technique was used in the present simulations. 
More elaborated methods may be more efficient, so that 
our results can be considered as lower limits on the ac¬ 
curacy that can be actually achieved using a given set of 
light curves. 

Although the Pelt method is well known, we briefly re¬ 
view the main steps followed to determine the time delay. 

We first define an interval of time delays 
[Atmin, Atmax], which Contains the true value of 
the time delay Atm. Note that with real data, predicted 
time delays for lenses are accurate enough to follow this 
approach, especially in cases where the redshifts of the 
lens and source are known. We then dehne N^ equally 
spaced time delays over the range [Atmin, Atmax], with 
interval < O.f days i.e. N^ > {Atmax - At™„)/0.1 -I- 1. 
The interval is small compared with the input time delay 
Atin and ensures that the precision of the results, even 
for Atin = 5 days, is not limited by the time resolution 
adopted in the simulations. 

The light curve of image B is then shifted, successively, 
through the set of Nd time delays. At. The problem is to 
find which curve B(At) best matches curve A, within the 
overlap region. For any curve B(At) the overlap region is 
defined as the set of points for which there exist points 
in curve A, both before and after in time. Curve A is 
then linearly interpolated to these points, and the disper¬ 
sion D‘^{At) in the magnitude differences between the two 
curves provides the measure of goodness of fit. Data points 
for which the distance from the interpolated date to the 
closest date in curve A, is greater than some specified limit 


(i.e. where the interpolation is unreliable) are ignored in 
this calculation. The search is limited, obviously, to time 
delays for which there is overlap of the two curves. Time 
delays of the order of half a year are thus only accessible 
for circumpolar objects. 

This procedure yields the dispersion spectrum D^{At). 
The true time delay Atout between the images should be 
evident as a minimum in the dispersion spectrum D^{At). 
This minimum is determined by htting a parabola to the 
dispersion spectrum. 

5. Results 

For every time delay Atin to be simulated, we explored the 
full range of 36 different combinations of the three param¬ 
eters, detailed in section 3. For each combination, we ran 
100,000 simulations, each time changing the quasar light 
curve, and modifying the observed points by adding ran¬ 
domly distributed errors (i.e. normally distributed devi¬ 
ates of 0.01 mag variance). The results for Atin = 80 days 
are summarized in Fig.|21 where the 36 panels correspond 
to the different parameter combinations. In each panel 
the measured time delays of the simulations are plotted 
in histogram form, with the measured mean and standard 
deviation (established by computing the range containing 
68 % of the results) quoted. The histograms are mostly 
symmetrical about their mean value, indicating that no 
strong systematic error is introduced. The slight shift (0.5 
days in the worst case) of the mean of the histogram rela¬ 
tive to value of Atm is small compared with the width of 
the histogram, i.e., the random error dominate the error 
budget. 

The results for the percentage error, for the 36 param¬ 
eter combinations, for all the simulated time delays, are 
presented in a compact way in Fig. ^ Here each panel 
shows the results for 12 parameter combinations, which 
are the 4 temporal samplings x the 3 visibilities, and the 
three columns correspond to the 3 amplitudes. Each row 
is for a different time delay. In this figure are also shown 
the values of the quoted relative errors for six published 
time delays of quasars monitored in optical wavelengths. 
These values are summarized in Table ^ and have been 
plotted at a sampling of 15 days, which roughly corre¬ 
spond to their effective sampling. Although the observa¬ 
tional strategy used for these quasars is not strictly the 
same as any of the ones we have defined, the predicted 
errors lie rather close to the real ones, with the excep¬ 
tion of the double quasar HE 2149-2745, but this quasar 
had very smooth variations over the two years of observa¬ 
tions, much smoother than the typical variations used in 
our simulations (see Fig. El). We do not include the twin 
quasar Q 0957-1-561 because it has a 10-year long light 
curve, much longer than the two years considered in our 
simulations. Moreover its time delay of 423 ± 6 days (Pelt 
et al. mm is much larger than the highest time delay 
used (i.e. 120 days). 

Although the predicted relative errors on the time de¬ 
lay are very close to the published values, they are usually 
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Atjjj = 80 days (2 years of observations) 
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Fig. 3. Histograms exploring the observational parameter space described in the text for the determination of a time 
delay of 80 days. Each curve is the probability density function for the time delay, obtained from 100,000 simulations, 
for a particular combination of the three variables. These are: 1- sampling interval, four columns, from left to right: 
irregular, 15 days, 7 days, 3 days); 2- visibility period, three bands from top to bottom: 12, 8, and 5 months; 3- peak-to- 
peak variation. A, three rows within each band, from top to bottom,: 0.3, 0.2, 0.1 mag. Each panel is labeled with the 
mean and standard deviation of the measured time delay, as well as the percentage error. The effect of microlensing 
is not included in these simulations, and is treated later. 


slightly more optimistic than the measured relative errors. 
The small discrepancies can be explained by differences in 
the parameters we use, compared with the characteristics 
of actual monitoring data, e.g.: 


— a shorter or longer monitoring period than the sup¬ 
posed two years, 

— different photometric errors than the supposed 0.01 
mag, 

— a temporal sampling that differs in detail from our 
idealized scheme. 


The most likely explanation for the simulations being 
too optimistic however remains the presence of microlens¬ 
ing in the light curves of real quasars, which is the subject 
of the next section. 


6. Influence of “slow” microlensing 

Not all the photometric variations observed in the light 
curves of the quasar images are intrinsic to the quasar. 
Microlensing by stars in the lensing galaxy introduces am¬ 
plification events that contaminate the light curves. 
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Fig. 4. Summary of the estimated percentage error on the measured time delay as a function of the observational 
parameters: 1- peak-to-peak variation, A; 2- sampling interval (x-axis of each panel); 3- visibility period. Each panel 
corresponds to one value of the input time delay The percentage error on the time delay, plotted on the y-axis, is 
calculated from 100,000 simulations. The lines connecting the points correspond to different periods of visibility. The 
solid lines are for circumpolar objects, the dotted lines are for the 8-month visibility, and the long dashed lines are 
for the 5-month visibility. We have used three peak-to-peak values of the amplitudes A for the simulated light curves, 
which increases from left to right in the three columns. The curve for the 5-month visibility and 120-day time delay 
has not been computed, as there are almost no data points in common between the curves of image A and image B. 
The star-shaped symbols plot the percentage errors for quasars with real measured optical time delays published in 
the literature. See text for further details. 


The severity of such events depends not only on the 
location of the images relative to the lens but also on 
whether the image considered is a minimum, maximum 
or a saddle point in the arrival time surface (Schechter & 
Wambsganss l2()()2'll . Consequently, the image closest to the 
lens, in projection on the plane of the sky, and hence with 
the larger density of potential microlenses, is not neces¬ 
sarily the one with more microlensing. The doubly lensed 
quasar HE 1104-1805 is a typical example, where the im¬ 


age the further away from the lens is the one with the 
largest microlensing events. 

Microlensing can act on different time scales, “slow” 
or “fast”, as compared with the time scale of the quasar’s 
intrinsic variations. A nice example of fast microlensing 
has been found in the light curve of HE 1104-1805 (e.g., 
Schechter et al. 1200311 . Since the temporal sampling used 
in past quasar monitoring programmes has been rather 
sparse, there is no other known example of fast microlens- 
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At;^ = 80 days (2 years of observations) 
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Fig. 5. Histograms exploring the observational parameter space described in the text for the determination of a time 
delay of 80 days, including the effects of microlensing. Each curve is the probability density function for the time 
delay, obtained from 100,000 simulations, for a sampling interval of 3 days, and for a particular combination of the 
variables. These are: 1- microlensing amplitude, = a ■ A, four columns, from left to right: a = 0,0.01,0.05,0.1; 
2- visibility period, three bands from top to bottom: 12, 8, and 5 months); 3- peak-to-peak variation, A, three rows 
within each band, from top to bottom: 0.3, 0.2, 0.1 mag. Each panel is labeled with the mean and standard deviation 
of the measured time delay, as well as the percentage error. While no strong systematic drifts of the histograms are 
seen relative to the input time delay Atm = 80 days, the width of the histograms are significantly broadened as 
microlensing increases. 


ing. Slow microlensing, with smooth variations spanning 
several months or even years are more common, or are at 
least better detected with existing data. The slow vari¬ 
ations in the four images of the Einstein Cross are the 
clearest examples of slow microlensing (e.g., Wozniak et 
al. l2000ll . 

Since most quasars with measured time delays have 
been shown to be affected by slow microlensing, it is 
mandatory to introduce this effect into our artificial light 
curves and to estimate how the time-delay measurement 


is modified. The slow microlensing events can be simu¬ 
lated by creating a microlensing light curve in the same 
manner as we did for the intrinsic variations of the quasar 
(i.e. by using a smooth random walk model), but with a 
different length scale and amplitude. We express the peak- 
to-peak microlensing amplitude A^ as a simple function 
of the quasar amplitude. We take it as A^ = a ■ A, with 
a = 0.01, 0.05, 0.10, in order to mimic a microlensing am¬ 
plitude of respectively 1, 5 and 10% of the amplitude of the 
quasar light curve. The microlensing curve is smoothed us- 
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Table 1. Published time delays and l-cr uncertainties 
for four lensed quasars, measured from their optical light 
curves. The percentage errors are given in parentheses. 


Object 


Time delay [days] 


Reference 


RXJ 0911+0551 
FBQ 0951+2635 
PG 1115+080 
SBS 1520+530 
B 1600+434 
HE 2149-2745 


146 ±4 (2.7%) 
16 ±2 (13%) 
23.7 ±3.4 (14%) 
130 ±3 (2.3%) 
51 ±2 (3.9%) 
103 ±12 (12%) 


Hjorth et al. I200(SI 
Jakobsson et al. 
Schechter et al. mm 
Burud et al. 12002^ 
Burud et al. 
Burud et al. 120021^ 


ing a convolution kernel of 100 days, i.e. ~ 3 times slower 
than the intrinsic variations of the quasar, and is added to 
the light curve of one of the quasar images. The choice of 
this image is irrelevant, because only relative differences 
between the lightcurves of the two components are con¬ 
sidered to extract the time delay. 

We then repeat the cross-correlation analysis. The 
microlensing event, thus, acts as an additional source 
of noise. Fast microlensing is not considered here. 
Introducing it is equivalent to artificially increase the 0.01 
mag error bar on the individual points. 

Fig. El plots the results for the case of the 80-day time 
delay, and 3-day sampling, with different amplitudes of mi¬ 
crolensing. The format is the same as in Fig. El It can be 
seen that no strong systematic variations are introduced. 
In each case the returned time delay is consistent with 
the input value, but microlensing substantially increases 
the uncertainty in the measurement, i.e. broadens the his¬ 
tograms. No distortion, i.e., skewness is apparent. The er¬ 
ror on the time delay measurement without microlensing 
(left column) typically degrades by a factor of approxi¬ 
mately 2 when 5% microlensing is present (i.e., a = 0.05), 
and by up to a factor of 6 with 10% microlensing. However, 
the shift between the mean of the distribution and 
is not much larger than in the no-microlensing case. Slow 
microlensing does not seem to introduce significant sys¬ 
tematic errors. 

Figs. El and 0 summarize all the results of our simu¬ 
lations including microlensing, in a way similar to Fig. 0 
showing how the error on the time delay degrades with in¬ 
creasing microlensing amplitude (a, plotted on the x-axis 
in each of the column plots). The figures are constructed 
for the irregular sampling as well as for the regular 3-day, 
7-day and 15-day samplings. 

We note that the time-delay determination is much 
more affected by microlensing with the 3-day sampling 
than with the 15-day or the ’’irregular” samplings: while 
changing a from 0 to 0.1 increases the relative error by 
a factor of 8 for the 3-day sampling, it increases only by 
a factor of 3-4 with the 15-day or the irregular sampling. 
Microlensing has a stronger effect on well sampled light 
curves than on sparser samplings. 


Similarly a light curve with large amplitude A will see 
its accuracy on the time-delay measurement slightly more 
degraded than one with a smaller amplitude. In both cases 
this may simply be due to the enhanced signal-to-noise of 
the light curves, either because more data points are avail¬ 
able, or because the quasar variations are stronger with 
respect to the photon noise of the individual photometric 
points. 

In general we can conclude that the more accurate the 
time-delay determination is in the case without microlens¬ 
ing, the more it degrades when a given amount a of mi¬ 
crolensing is added: better data are more sensitive to mi¬ 
crolensing. On the other hand, the data allowing accurate 
time delay determinations in the absence of microlensing 
are usually also better suited to the accurate subtraction 
of the microlensing events. 

7. When the time delay becomes close to the 

length of the visibility window 

So far we have compared light curves sampled with regu¬ 
lar samplings, plus one irregular sampling. The main dif¬ 
ference between these samplings was the number of data 
points within the period of 2 years of observations. It is 
then not surprising that finer sampling leads to better re¬ 
sults. The simulations we have done allow us to quantify 
the error bar on the time delay for each sampling. 

Another natural question arising is: is there an opti¬ 
mal way to distribute a fixed number of sampling points, 
in order to reach the best possible accuracy on the time 
delay ? This has been explored in other areas of astron¬ 
omy, for example by adopting a logarithmic sampling of 
the data. We have tested the effect of such a sampling 
on quasar light curves. Fig. |H1 shows the results of the 
simulations, where we compare the (regular) 7-day sam¬ 
pling to a sampling adopting the exact same number of 
data points but distributed in a logarithmic way. As for 
the regular case, we have introduced a small randomly 
distributed error (±0.4 days) on each observing data to 
account for weather conditions and scheduling. As shown 
in Fig. 0 the curve starts with a very high frequency of 
observations and continues with a sampling getting close 
to regular. An important consequence is that objects that 
have a time delay of the order of the visibility period will 
be well sampled exactly where the two quasar light curves 
significantly overlap after correcting for the time delay. In 
other words, the logarithmic scale allows to sample very 
well the (short) parts of the curves that will overlap after 
the time delay is applied. 

The result in Fig. 0 is striking. As soon as the time 
delay is close to the length of the visibility window, the 
regular method fails to produce symmetrical histograms, 
whereas the histograms obtained with the logarithmic 
scale are narrower and more symmetrical about their 
mean. Their mean is also closer to than with the reg¬ 
ular method. This is no longer true when the time delay is 
shorter than the visibility window, where the logarithmic 
sampling even degrades the results. 
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Fig. 8. Time delay distributions for four different values of each time for 3 visibilities and a peak-to-peak 

variation A = 0.3. We compare the distributions obtained for the logarithmic sampling (without microlensing), with 
the results for the 7-day sampling. Clear distortions of the histograms are seen when the time delay is close to the 
visibility window of the object, when a regular sampling is adopted (left column in the four panels). The histograms 
obtained with the logarithmic sampling are well symmetrical and narrow. The systematic error is also reduced. This 
effect is not so evident when the time delay is much shorter than the visibility window, where the logarithmic sampling 
even degrades the results. 


8. Conclusions 

We have undertaken a set of simple but realistic numerical 
simulations in order to optimize the observing strategy 
of our COSMOGRAIL photometric monitoring programs 
aimed at measuring Hq. The predicted error bars on time 
delays compare very well with the ones obtained in optical 
wavelengths with real data. 

It is immediately seen from Figs. IHl and [7| that short 
time delays will never be measured accurately, i.e., with a 
precision better than 2%, unless the quasar amplitude A is 
substantially larger than 0.3 mag. Even with no microlens¬ 
ing and the 3-day sampling, time delays shorter than 10 
days are measurable with 10% accuracy, at best. Time 
delays between 40 and 100 days seem optimal, especially 


in the case of circumpolar objects, where one can easily 
achieve 2% accuracy, even with the 7-day sampling. 

Equatorial objects should be avoided. Although they 
are accessible from the north and south, they are visible 
under good conditions for only 5-6 months along the year. 
This makes it impossible to measure time delays larger 
than 100 days (hence the corresponding long-dashed curve 
is not represented in the relevant figures). For shorter time 
delays, e.g., 80 days, the estimated error for an equatorial 
object is twice that of the same object if it were circum¬ 
polar. 

Microlensing complicates the situation. With 5% mi¬ 
crolensing (as defined here), the estimated error on the 
time delay is twice that of the no-microlensing case. Again, 
optimal time delays are around 100 days, with a visibil- 
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ity of at least 8 months. Assuming an amplitude A = 0.2 
magnitude and 5% microlensing, an accuracy of 2% on 
the time delay is still possible for these objects. The long 
time delays also allow a sampling as long as 7 days to be 
adopted. 

While microlensing increases the random error on the 
time delay, it does not increase significantly the systematic 
error (i.e., \Atin — Atoutl), which remains low, usually 5 
to 10 times lower than the random error, with or without 
microlensing. 

Finally, adopting a logarithmic sampling step can im¬ 
prove the accuracy on the time delay when the time delay 
is close to the length of the visibility window of the ob¬ 
ject. However, this logarithmic sampling usually degrades 
the results for all other combinations of time delays and 
visibility windows. 

The game of defining what could be a “golden” lens, at 
least in terms of quality of the time delay measurement, 
is not an easy one. This is why we have attempted in 
this paper to provide a grid of predicted error bars on the 
time delay, based on simple assumptions. The results are 
presented in a compact way in Figs. El and Figs. 0 We 
will use these plots to choose optimal combinations of 1 - 
predicted time delay, 2- object visibility and 3- temporal 
sampling, to reach a target accuracy on the time delay. 
Even with large amounts of telescope time, the number of 
new lensed quasars is increasing quickly and one has to 
select the best possible cases. We hope that the present 
work will be useful for the task of identifying the objects 
that are the most likely to be measured accurately, so that 
the only significant remaining source of uncertainty on Hq 
will be the lens model. 
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